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The flavor asymmetry of the light quark sea of the nucleon is determined in the kinematic range 
0.02 < x < 0.3 and 1 GeV 2 < Q 2 < 10 GeV 2 , for the first time fr om semi-inclusive deep-inelastic 
scattering. The quantity (d(x) — u(x))/(u(x) — d(x)) is derived from a relationship between the 
yields of positive and negative pions from unpolarized hydrogen and deuterium targets. The flavor 
asymmetry d — u is found to be non-zero and x dependent, showing an excess of d over u quarks in 
the proton. 
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The flavor content of the nucleon sea has come to be 
recognized as an important domain for testing models of 
nucleon structure jlj] . Until recently there has been little 
experimental constraint on the flavor asymmetry in the 
quark distributions of the light sea quarks in the nucleon. 
The first was an integral test, based on a comparison of 
inclusive deep-inelastic scattering on the proton and the 
neutron to determine the Gottfried Sum H defined as 



S G 



1 dx 



(F£(x)-F?(x)), 



(1) 



where FJ" (x) and F%{x) are the structure functions of 
the proton and neutron, respectively. The assumption of 
isospin symmetry in the quark-parton model allows the 
Sum to be written as 



5, 



1 



G 



(u v (x) — d v (x))dx 



(d(x) — u(x))dx, 
(2) 



where u v (x) and d v (x) are the density functions of the 
valence quarks and u(x) and d(x) those of the antiquarks 
in the proton, and x is the fraction of the nucleon light 
cone momentum carried by the struck quark. A flavor 
symmetric sea, u(x) = d(x), leads to the Gottfried Sum 
Rule Sg = 1/3. A measurement of the Gottfried Sum by 
NMC § resulted in S G = 0.235±0.026, at Q 2 = 4 GeV 2 . 
If isospin symmetry holds, a global flavor asymmetry 

(d(x) — u(x))dx w 0.15 would account for the NMC re- 
sult. Many ideas such as Pauli blocking and pion clouds 
have been proposed to explain such a flavor asymme- 
try in the sea Two methods have been proposed to 
measure its x dependence: the Drell-Yan process B and 
semi- inclusive deep inelastic scattering Q . 

Results are presented here for the x dependence of 
(d—u)/(u — d) from the analysis of charged pion yields 
in semi-inclusive deep-inelastic scattering of positrons 
on unpolarized hydrogen and deuterium targets. The 
data cover the kinematic range 0.02 < x < 0.3 and 



1 GeV 2 < Q 2 < 10 GeV 2 , where -Q 2 is the four- 
momentum transfer squared of the exchanged virtual 
photon and x = Q 2 / '(2 Mi/) is the Bjorken variable with 
M being the proton mass, v = E — E' the virtual photon 
energy, and E(E') the energy of the incident (scattered) 
positron. 

The ratio (d—u)/(u — d) is determined from a quan- 
tity that combines good sensitivity to the sea asymmetry 
with minimal sensitivity to instrumental effects. This is 
the ratio of the differences between charged pion yields 
for proton and neutron targets B: 



r(x, z) 



Nf(x,z)-NZ~(x,z) 
N^{x,z)-N- + {x,zY 



(3) 



where z = E 71 jv is the fraction of the virtual photon 
energy carried by the pion and N 7T (x,z) is the yield of 
pions coming from deep-inelastic scattering off nucleons. 
Factorization between the hard scattering process and 
the hadronization of the struck quark implies 



: (x, z)«^ej [qi(x)D^ (z) + q^D^ (z)] , (4) 



where is the quark charge in units of the elementary 
charge, qi(x) and ~q~i(x) are the density distributions of 
quarks and antiquarks of flavor i, and the fragmenta- 
tion functions D*. (z) represent the probability that the 
quark of flavor i fragments to a charged pion. Assum- 
ing isospin symmetry between protons and neutrons as 
well as charge conjugation invariance, the number of light 
quark fragmentation functions is reduced to two: the 
favored and disfavored fragmentation functions ||, e.g. 

and . Using Eq. [|to express r(x, z) in terms of 
quark distributions and fragmentation functions results 
in 



1 + r(x, z) u(x) — d(x) + u(x) — d{x) 



l — r(x,z) (u(x) - u(x)) — (d(x) — d(x)) 



J(z), (5) 
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where J(z) = f (0^) and D\z) = (z)/Dl + (z). 
It should be noted that the right-hand side of Eq. || fac- 
torizes into two independent functions of x and z respec- 
tively. Thus, the equation may be rearranged to isolate 
a quantity sensitive to the flavor asymmetry: 

d{x) - u(x) = J(z)[l-r(x,z)] - [l + r(x,z)} 
u(x) — d(x) J{z)[l — r(x, z)] + [1 + r(x, z)] 

The assumptions leading to Eq. imply that the expres- 
sion on the right side of Eq. ^ is independent of z. 

In the HERMES experiment, 27.5 GeV positrons cir- 
culating in the HERA storage ring at DESY are scat- 
tered on internal hydrogen (H2), deuterium (D2), and 
helium-3 ( 3 He) targets. The target gas is fed into the cen- 
ter of a 40-cm-long open-ended storage cell. Systematic 
uncertainties are minimized by cycling the three gases 
through one of various permutations, several times dur- 
ing each 8-12 hour HERA fill. The unpolarized lumi- 
nosity of (1 — 5) x 10 32 cm~ 2 s _1 is monitored by a pair 
of NaBi(WC>4)2 electromagnetic calorimeters that detect 
Bhabha scattering from target electrons. 

A detailed description of the HERMES spectrometer is 
provided elsewhere . It is a forward spectrometer that 
identifies the scattered positron as well as the associated 
hadrons, in the scattering angle range of 0.04 rad < 6 < 
0.22 rad. Positron-hadron discrimination is based on in- 
formation from four particle identification detectors: a 
threshold gas Cerenkov counter, a transition-radiation 
detector, a lead-glass electromagnetic calorimeter, and a 
preshower detector located directly before the calorime- 
ter. This combination provides an average positron iden- 
tification efficiency of 99%, with a hadron contamination 
that is dependent on the kinematics of the positron but 
is always less than 1%. One of the features of the spec- 
trometer central to this work is its threshold Cerenkov 
counter, which distinguishes pions from heavier hadrons 
for particles with momenta greater than 3.8 GeV while 
maintaining better than 95% pion identification efficiency 
at momenta larger than 6 GeV. 

The kinematic requirements used in this analysis are 
Q 2 > 1 GeV 2 , invariant mass W of the initial photon- 
nucleon system greater than 2 GeV, and fractional en- 
ergy y = v I 'E of the incident lepton transferred to the 
virtual photon less than 0.85. Events are required to 
originate from the target-beam interaction region. For 
the determination of r(x, z), the data are partitioned 
in five bins of x in the range 0.02 < x < 0.3 and up 
to six bins in z in the range 0.2 < z < 0.8. The pi- 
ons are selected to be in the current fragmentation re- 
gion (i.e. derived from the struck quark and not tar- 
get fragments) by requiring xp — 2Pp/W to be greater 
than 0.25, where Pj, is the momentum component of the 
pion in the longitudinal direction with respect to the vir- 
tual photon in the photon-nucleon center-of-mass frame. 
Data at Bjorken x > 0.3 are excluded because they are 



dominated by scattering from the valence quarks. The 
yields for the neutron target appearing in Eq. |^ can be 
inferred from yields obtained from deuterium and hydro- 
gen targets, N% = N% — N£ . Thus the determination 
of (d(x) — u(x))/(u(x) — d(x)) is made by taking ratios of 
differences between yields from hydrogen and deuterium. 
The pion yields extracted from the data are N% et {x, z) / C 
where NJ et (%, z) is the number of detected pions associ- 
ated with deep-inelastic scattering events, and C is the in- 
tegrated luminosity (nuclei/cm 2 /s) corrected for instru- 
mental dead-time and inefficiencies; the possible effect of 
different acceptances for 7r + and 7r~ is included in the 
systematic uncertainties discussed below. The fragmen- 
tation functions (z) and their ratio D'(z) were ex- 
tracted from the measured pion yields from each of the 
unpolarised targets: H 2 , D 2 and 3 He ||, using Eq. 4 
together with the GRV 94 LO parameterization of 
parton distributions. No statistically significant differ- 
ences were found between the results from the different 
targets. To avoid statistical correlations with the pion 
ratios, only D'(z) extracted from 3 He is used in this anal- 
ysis. The fragmentation function ratio D'(z) determined 
within the HERMES acceptance is found to vary with x 
and Q 2 typically by less than 10%, except in the low-yield 
region near z = 0.8 where as much as 40% may be seen. 
However, extraction of the flavor asymmetry in bins of x 
and Q 2 revealed the averaged result to be unaffected by 
this variation. The dependence of the left-right instru- 
mental asymmetry on pion charge is significant only at 
the lowest momenta and is taken into account as a cor- 
rection. For all x bins, the resulting change in the sea 
asymmetry is less than 1.5%. 

The results for the quantity (d — u) / (u — d) are 
shown in Fig. ^ as a function of z for five bins in x. 
The data show no z dependence and are consistent with 
the form of factorization shown in Eq. |[ The values of 
(d — u)/(u — d) averaged over z are plotted in Fig. ||(a) 
as a function of x, and are presented in Table Q. They 
are non-zero and positive everywhere in the measured 
x region, clearly showing an excess of d quarks over u 
quarks in the proton. Also included in this figure are the 
GRV 94 LO §, CTEQ 41q ||, MRS {A)_ jljj low Q 2 , 
and MRST (98) jl3| parameterizations of (d-u)/(u-d), 
calculated at the average Q 2 of the events in each x bin. 
The HERMES results are determined without making 
higher order QCD corrections and so only leading or- 
der parameterizations should be compared to them. At 
low Q 2 , only GRV 94 is available at leading order. The 
distribution d — u as a function of x is thus derived 
from (d — u)/(u — d) using the GRV 94 LO param- 
eterization of u(x) — d(x), again calculated at the ap- 
propriate Q 2 for each x bin. Other combinations from 
GRV 94 LO could have been used instead, such aS Zi V &y 
or u + u — d — d, but these choices result in more sensi- 
tivity to this model input, as indicated by differences 
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seen from substituting the GRV 98 LO parameteriza- 
tion The results are presented in Fig. ||(b) and Ta- 
ble JL and give an integral over the measured x region of 
fow(d(x) - u{x))dx = 0.107 ± 0.021(stat) ± 0.017(sys) , 
after the integrand at each x bin is 'evolved' to the av- 
erage measured Q 2 of 2.3 GeV 2 by applying the ratio of 
the GRV 94 values for d — u at the two values of Q 2 . 
The GRV 94 LO parameterization is used to estimate the 
portion of the integral in the unmeasured regions. The 
contribution from x < 0.02 is significant: 0.050 (almost 
50% of that over the measured x range), while the high 
x region x > 0.3 contributes only 0.006. Including an 
extrapolation uncertainty estimate based on differences 
among the available parameterizations, the total integral 
over all x is then approximately 0.16 ± 0.03, which is 
consistent with the results first seen by NMC j|. 

An alternate approach to the measurement of the fla- 
vor asymmetry is the comparison of the Drell-Yan process 
on protons and deuterons, which is sensitive to the quan- 
tity u/d Using this approach, the NA51 collaboration 
at CERN determined that u/d = 0.51±0.04±0.05 at x = 
0.18 and at a center-of-mass energy of a/s = 29 GeV p4| . 
The E866 collaboration at Fermilab recently reported 
a measurement of the flavor asymmetry in the x range 
0.02 < x < 0.345 ||_. E866 usesjthe CTEQ 4M pa- 
rameterization of u+d to extract d—u from the measured 
quantity d/u. These results given at a Q 2 of 54 GeV 2 are 
included in Fig. ||(b). This figure demonstrates that the 
sea asymmetry measured in deep-inelastic scattering and 
in Drell-Yan experiments agree, even though the Q 2 of 
the two experiments differ by a factor of about 20. 

The total systematic uncertainties associated with the 
determination of the quantity (d — u)/(u — d) in the 
present work are given in Table [|. These uncertainties 
include various effects as follows. After the correction at 
low momenta, any residual charge-dependent instrumen- 
tal acceptance asymmetry is estimated to be less than 
1.5%. Effects of the finite resolution of the spectrometer 
and bremsstrahlung have been investigated with Monte 
Carlo studies. These effects are small and are included in 
the systematic uncertainty. Radiative corrections to the 
pion yields have been calculated following the prescip- 
tion of reference Jr|. Their effect on the pion ratio in 
Eq. H is estimated to be less than 1% and is included in 
the systematic uncertainty. Lepton contamination in the 
pion sample is found to be at most of order one percent, 
and the differences between the two targets were no more 
than 0.2%. The effects of hadrons misidentified as pions 
have been estimated from the Monte Carlo simulations 
and are also included. Finally, the uncertainty in the 
fragmentation function ratio D'(z) as measured at HER- 
MES has been accounted for. The statistical uncertainty 
of D'(z) has been included in quadrature to the statisti- 
cal uncertainty of the pion yields. The systematic com- 
ponent of this ratio determined from systematic studies 



of the HERMES data has been included in quadrature 
to the total uncertainty. 

The only additional systematic uncertainty associated 
with d—u is the choice of parameterization to extract 
d-u from (d - u)/(u - d). The GRV 94 LO § pa- 
rameterization of u(x) — d(x) is used. The MRS (A) 0] 
low Q 2 parameterization is used as an alternative. The 
resulting difference in the value of d — u is used as an 
estimate of this uncertainty. 

In summary, the flavor asymmetry of the light quark 
sea is extracted for the first time from pion yields in semi- 
inclusive deep-inelastic scattering. The data are recorded 
at the HERMES experiment in the kinematic range of 
0.02 < x < 0.3 and 1 GeV 2 < Q 2 < 10 GeV 2 . The mea- 
sured quantity (d—u) / (u—d) is found to be non-zero over 
the entire x range measured, clearly showing an excess 
of d quarks over u quarks in the proton sea. Its x de- 
pendence will serve as a further constraint on global pa- 
rameterizations of parton density functions. The integral 
over the measured region of the derived quantity d—u is 
J°^ 2 (d(x)-u{x))dx= 0.107±0.021(stat)±0.017(sys) at 
a Q 2 of 2.3 GeV 2 , which accounts for two-thirds of the 
Gottfried Sum Rule deficit. 
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TABLE I. Values for (d - u)/{u - d) and d-u. 
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range 


< x > 


<Q 2 > 
(GeV 2 ) 


d — u 
u — d 


Stat, 
error 


Sys. 
error 


d — u 


Stat, 
error 


Sys. 
error 


0.020-0.050 


0.038 


1.33 


0.20 


0.08 


0.07 


0.53 


0.21 


0.20 
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0.062 


1.82 


0.42 
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0.91 
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0.15 
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0.131 
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0.08 


0.05 


0.36 


0.13 


0.10 
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0.198 


4.88 


0.17 
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0.11 


0.09 
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FIG. 1. The distribution (d — u)/(u — d) as a function of z in five bins of x. The points are fit to a constant for each x bin. 
The error bars represent statistical and systematic uncertainties added in quadrature. Some points are omitted due to limited 
statistics. 
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(u -d) as a function of x. Also included are the GRV 94 LO §, CTEQ 41q [0, MRS (A) Q low 
parameterizations calculated at the appropriate Q 2 for each a;- bin. (b) d — u as a function of a;. The 



FIG. 2. (a) (d - u 
Q 2 , and MRST (98) _ 

open circles represent the E866 Q determination of d — u at Qo=54 GeV 2 . The solid (dashed) curve is the GRV 94 LO ]S|] 
parameterization, evaluated at the Q 2 values for the HERMES (E866) experiment. The inner error bars represent the statistical 
uncertainties while for all data including those from E866, the total error bars represent statistical and systematic uncertainties 
added in quadrature. 
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